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What is attitude? 


• ‘Attitude' means the orientation of a coordinate 
Irame fixed in the spacecraft body relative to a 
reference coordinate frame. 




What do we observe? 

Angles or vectors to known references 

Sun (coarse sun sensor or digital sun sensor) 
harth (horizon sensor or landmarks) 

Stars (star trackers or fine guidance sensors) 
Magnetic field (triaxial magnetometers) 

— RF sources, like GPS (phase interferometry) 

- Motion (inertial sensors — gyros) 








Euler's Theorem (1775) 

The general displacement 
of a rigid body with one 
point fixed is a rotation 
about some axis 
r in reference frame, 
b in body frame 
b = 7?(e, 0) r, where 
R(e , 0) = ee T + cos 0(7- ee T ) - sin0[ex] 
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Cross Product Matrix 

• [ex] is the cross product matrix 


[ex]s 


Defined so that [ex]v = exv 





Rotation Matrix 


• The 3x3 rotation matrix R is 
-orthogonal (transpose - inverse), and 

proper (determinant = f 1 ) 

• Also called the direction cosine matrix (DCM) 

• Also called the attitude matrix, denoted by A or D 

• Attitude kinematics equation dA/dt = -[©x ]A , 
where G) is the angular velocity or body rate vector 

• Specification of attitude requires three parameters 




Euler Angles 

• Parameterization as product of three rotations 
R ljk {& 9, y/)=R(t k , y) R(e r 0) R(e n 0) 

- Symmetric, ijk — 131, 121,232,212,313,323 
-Asymmetric, ijk — 123, 132,231,213,312,321 

- ‘Gimbal-lock 1 singularity for some 0 

• Generalized Euler Angles (Davenport, 1973) 

- Three rotation axes are not coordinate axes 

- Must have tj perpendicular to e i and e A 


'NASA, 

Other Parameterizations 

• It is topologically impossible to have a global 
nonsingular 3-dimensional parameterization 

• Singular 3-dimensional parameterizations 

- Rodrigues parameters, also known as Cayley 
parameters or the Gibbs vector = etan(0/2), 

- Modified Rodrigues parameters, etan(0/4) 

• Nonsingular 4-dimensional parameterization 
-Euler-Rodrigues parameters, or quaternion 

• ,..****. if . ' 





Olinde Rodrigues (1795-1850) 


• Son of a Jewish accountant in Bordeaux 

• Doctorate in math in 1816 from U. of Paris 

- Rodrigues formula for Legendre polynomials 

• Banking and utopian socialism for 24 years 

Important in introducing railroads to France 
Edited an anthology of workers’ poetry 

• Published his seminal paper on analysis of 
rotations in Liouville’s Journal in 1840 



Why Quaternions? 

• 4-dimensional non-singular parameterization 
q = [</, q 2 t h f e si n(0/2). q 4 - cos(0/2) 

One constraint. |c/| 2 = 1 

• Rotation matrix, using half-angle formulas, is 

R (<? ) ’ Uh 2 ~ |q| 2 ) / + 2qq T - 2 q 4 [q x] 


• Simple product rule R(q)R{q') = R{q®q') 

• Quaternion kinematics d q/6t = h (D ® q 
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attitude Determination 
Methods 

‘Single Frame’ Methods 

- Ad hoc methods 

- TRIAD (Black, 1964) 

-Optimal methods (Wahba’s problem, 1965) 
Filters 

- Extended Kalman Filter (EKF) 

- Hqq filter (Berman, Markley & Shaked, 1993) 

- Predictive filter (Crassidis & Markley, 1997) 




Hubble Space Telescope (HST) 
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Launch April 1990 
Mass = 1 1000 kg 
Pointing = 0.007 arc seconds 
Attitude sensors: 

- Fine Guidance Sensors 

- Gyros 

- Star trackers 

- Sun sensors 

- Magnetometers 






HST Science 
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Solar, Anomalous, and 
Magnetospheric Particle Explorer 
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• Launch July 1992 

• The first Small Explorer 
(SMEX) mission 

• Mass = 160 kg 


Pointing = 2 ‘~ 


Attitude sensors 

- Sun Sensors (0.25°) 

- Magnetometers 




Small Explorer Satellite Finds Radiation Belt 
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Two Contrasting Missions 

Differences 

- HST pointing six orders of magnitude finer 

- HST more complex, uses more sensors 
Similarities 

- Onboard attitude determination 

- Programmable digital computer 

- HST gyroless safemode uses the same sensor 
information as SAMPEX 




SAMPEX Attitude 
Determination 
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• TRIAD is used for attitude determination 
-A ‘single frame’ method 

- Measurements at one instant of time 

• Measurements represented as (two) unit vectors 

- Direction from the spacecraft to the Sun 

- Direction of the Earth’s magnetic field 

- Coordinates in reference and body frames 

- Minimum number required for a solution 



TRIAD (Black, 1964) 

• Orthogonal triad of body frame unit vectors: 
b, - spacecraft-to-Sun unit vector 

b 2 = perpendicular to Sun and magnetic field 
b 3 = b, x b 2 

• Corresponding triad of reference frame vectors 
r,, r 2 , r 3 from ephemerides and field models 

• A = [ b, b 2 b 3 ][r, r 2 r 3 f - 1, b, r, T 





Wahba's Problem (1965) 


• Optimal estimate from n vector observations 

• Find the proper orthogonal A that minimizes the 
loss function, with positive weights a n 

m-teafo-Arp, 

where {b,} are the body frame unit vectors, and 
{r,} are the same vectors in the reference frame. 

• Can write L(A) = traced# 1 ), where 

B = Z,a,b, r, T 
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Procrustes Problem 

Wahba’s problem is equivalent to the Procrustes 
problem: to find the orthogonal matrix A that is 
closest to B in the Frobenius norm, defined by 

|| M\ | 2 = M fJ 2 = trace(A/A/ T ) 

\\A -B\\ 2 = traced - B)(A -B) T ] 

= trace(^ T )-2trace(^ £ T )+trace(£5 T ) 
Since trace(A4 T ) =trace( 7) = 3, 
this is minimized by ma ximizing traced 7? T ) 





Davenport's q Method (1977) 

AUl)~Ul A - |q|-)/+2qq T 2r/,[qx]. 


is a homogeneous quadratic function off/, so 
UA) = I, a, - trac q(AB t ) - ] Z l a l - q T Kq. 

where K is a symmetric, traceless 4x4 matrix 

whose elements are linear in the elements of B. 

The optimal quaternion <y opt is the eigenvector of 
K with the maximum eigenvalue 


NASA 


Faster (but less robust) 
Solutions of Wahba's Problem 

Shuster s QlJatemion ESTimator (1978) 

- Solve the equation det (A / - K) = 0 iteratively, 

starting from I, a, = L(A opl ) 

- Then simple matrix algebra gives q 

- Error covariance is P = [ Z # a, (/ - b, b, T )] 1 

if weights are inverse measurement variances. 
Mortari’s ESOQ (1996) and ESOQ2 (1997) 
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for Solving 
Wahba’s Problem 

• Singular Value Decomposition Method (1987) 
B = U diag([.v h 5 2 , .v 3 ]) l r] \ then A opt = UV 7 
with U, V proper orthogonal and 5 , > s 2 > |s 3 1. 
P= U diag([(s 2 +s 3 ) ~ l Xs 3 +s l )~\(s l +s 2 ) ~'])£/ T 

• Fast Optimal Attitude Matrix (FOAM, 1992) 
^o P t = C 1 [(*+1 W )B + \ iax adjfl + BB T Bl 

with £= = and ^=i(^ iax 2 - M 2 ) 

P=£ 1 (kI+BB j ) 


My Algorithms 




Tropical Rainfall Measuring 
Mission (TRMM) 

• Launch November 1997 

• Mass = 3500 kg (largest 
jH?* ever built at GSF C) 

• Pointing = 0.2° 

• Attitude sensors: 

- Static Earth sensor 

- Gyros 

- Sun sensors 
ometers 
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TRMM Attitude Determination 

• Science mode uses R ilx {ya\\\ pitch. roll) 

- Relative to an Earth-pointing reference 

- Earth sensor gives pitch and roll directly 
-yaw determined from gyro and sun sensors 

- Gimbal lock would be at pitch = 90° 


• Independent contingency mode 

- Uses gyros, sun sensors, and magnetometer 

- Quaternion EKF 







(MAP), a MIDEX 



Launch November 2000 
L2 orbit, 1.5 million km 
behind the Earth 
Mass = 830 kg 
Pointing = 0.03° 

Attitude sensors: 

- Star trackers 

- Gyros 

- Sun sensors 
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Quaternion EKF 

• Straightforward "additive’ EKF has q = q est + A q 

• EKF uses error covariance P 4x4 = E{(A<y)(A<y) T } 

• The normalization constraint is, to first order, 

1 = M 2 = k/estl 2 + 2( Ary) T cy est = 1 + 2(Acj) T q cst . 

• So (A q) T q cs = 0 to first order in the error. 

• Thus P 4x4 q est = E { (Aq)(Aq) T ) <y est = 0. 

• The covariance is not positive definite. 

• What to do? 


Multiplicative Quaternion EKF 

q = 5cj ® g est , with Sej and q est unit quaternions. 

Use a local nonsingular 3-dimensional 
parameterization x of Sq ( e.g . the Gibbs vector). 

The multiplicative EKF filters x instead of q. 

P 2 x 3 - E{xx T } is nonsingular. 

After an update, q est is reset with x information. 
This method was invented in 1969 (SPARS), and 
has been used in NASA spacecraft since 1978. 
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Swift Gamma Ray Burst 




Explorer, another MIDEX 



• Launch 2003 

• Mass - 1 300 kg 

• Pointing - 3 arc seconds 

• Attitude sensors: 

- Star trackers 
Ring Laser Gyros 
Sun sensors 

- Magnetometer 





Summary 



Attitude iequires 3 parameters, but there is no 
global nonsingular 3-d parameterization 


• Onboard attitude determination is the norm, 
using either ‘single frame’ or filtering methods 

• Several mathematical representations of attitude 
are used, often on the same spacecraft; 
quaternions are used in the EKF and QUEST 


Attitude detennination precision requirements 
vary over at least six orders of magnitude 





